Spectroscopic Evidence for Anisotropic S-Wave Pairing Symmetry in MgB 2 
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Scanning tunneling spectroscopy of superconducting MgB2 (T c = 39 K) were studied on high- 
density pellets and c-axis oriented films. The sample surfaces were chemically etched to remove 
surface carbonates and hydroxides, and the data were compared with calculated spectra for all 
symmetry-allowed pairing channels. The pairing potential (Afc) is best described by an anisotropic 
s-wave pairing model, with Afc = A xy sin 2 6k + A z cos 2 #fc, where 6k is the angle relative to the 
crystalline c-axis, A z ~ 8.0 meV, and A xy ~ 5.0 meV. 



Since the discovery of superconductivity in MgPj2 at a 
superconducting transition temperature T c ~ 39 K [Q, 
a number of reports Jl|, ||, |], |], ||, ||, fjj have suggested 
that this hole-doped layered superconductor ||, ||, 
may be consistent with conventional BCS s-wave pairing. 
On the other hand, muon spin rotation (/iSR) studies of 
MgB 2 have found that the temperature dependence of 
the magnetic penetration depth is suggestive of uncon- 
ventional pairing symmetry with nodes in the supercon- 
ducting order parameter [jy]. To address the issue of 
the pairing symmetry in this new superconductor, possi- 
ble complications by disorder or surface impurities must 
be considered. Indeed, recent x-ray photoemission spec- 
troscopy (XPS) studies have revealed that MgC03 and 
Mg(OH) 2 exist on the surface of as-grown MgB 2 p2| . It 
is therefore important to understand how these surface 
impurity phases may contribute to surface-sensitive ex- 
periments such as the scanning tunneling spectroscopy 
(STS) [H H U and point-contact measurements of 
the quasiparticle spectra. In particular, existing STS 
data on as-grown polycrystalline MgB 2 || exhibited "V- 
shape" differential conductance (dl^s /dV) versus volt- 
age (V) plots near zero-bias (i.e. the Fermi level Ef), 
with rounded "humps" rather than sharp peaks at the 
gap values (V = ±A/e) and large residual density of 
states (DOS) at Ep. Those spectra were fitted with an 
s-wave pairing potential A broadened by disorder param- 
eterized as T, and a large ratio of (r/A) ~ 60% was sug- 
gested S. For comparison, in cuprate superconductors 
the V-shape conductance spectra near Ep for quasiparti- 
cle tunneling along the c-axis are known to be the signa- 
ture of the d x 2_ y 2 pairing symmetry [ll| 16, p"7|, 
and strong directionality in the quasiparticle spectra has 
been observed (ll| |l5l. I n particular, a zero-bias con- 
ductance peak (ZBCP) Jl8[ fl9| can occur if quasiparti- 
cles are incident close to the {110} nodal direction of the 
d2.2_ y 2-wave order parameter. Thus, should the pairing 
symmetry be unconventional, the observation of V-shape 
tunneling spectra in polycrystalline MgB 2 samples asso- 



ciated with certain grain orientations would be accom- 
panied by frequent occurrence of ZBCP for other grain 
orientations. To date, no ZBCP has been found from vac- 
uum tunneling studies of as-grown MgB 2 H, ||, ||] . How- 
ever, a major concern presented by existing quasiparticle 
spectra is that the measured gap values vary widely, and 
that most values are smaller than that the BCS predic- 
tion|,|,lJ]. 

Our starting point for investigating the pairing sym- 
metry of MgB 2 is to consider all the possible pairing 
channels based on group theory. The global symmetry 
group Q of MgB 2 in its normal state can be expressed by 
Q = U(l) x T x SU(2) x g space , where U(l) is the elec- 
tromagnetic gauge broken below T c , T and SU (2) denote 
the time-reversal and spin-rotational symmetries that 
are generally preserved below T c for spin-singlet Cooper 
pairs, and Gspace is the space group Dq^ for MgB 2 . Given 
that the Cooper pairs in MgB 2 are spin-singlets J3 and 
that no other obvious symmetry-breaking fields exist be- 
low T c except U(l), the possible pairing channels can 
be derived from the even-parity irreducible representa- 
tions of D$h- For a single-component superconductor, 
the relevant pairing channels can be further reduced to 
four one-dimensional (ID) even-parity irreducible repre- 
sentations in Dgh- Ai g , A 2ff , B\ g and -B 2g . The pairing 
potentials A^ for these representations can be expressed 
as a function of the momentum k to the lowest order: 
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is 



Afc = A , (isotropic s) 
Afc = A [l + ecos(60fc)], (anisotropic s) 
Afc = A xy sin 2 6k + A 2 cos 2 6k, (anisotropic s) 

A fe = A o sin 6 (9 fe sin(60 fe ), (1) 
A k = A cos 6 k sin 3 6 k sin(3<f> k ), 
Ak — A cos 6k sin 3 6k cos(3</>/c). 



Here 6k is the angle measured relative to k z , with k z 
parallel to the crystalline c-axis, and 4>k is measured rel- 
ative to k x . In addition, < e < 1 and A xy ^ A z for 
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FIG. 1: Right panels: Graphical representations for possible 
order parameters permitted by the group symmetry and 
spin-singlet pairing. Left panels: Simulated differential con- 
ductance (Gjvjvd/jvs /dV) vs. voltage (V) quasiparticle tun- 
neling spectra at 4.2 K, assuming Ao = 6.5 meV, for the fol- 
lowing ID even-parity representations, (a) Ai g , anisotropic 
s-wave with in-plane anisotropy; (b) Ai g , anisotropic s-wave 
with uniaxial symmetry; (c) A2 g ; (d) Bi g , or Bi g by rotating 
Bi g order parameter through an angle (it/ 6) relative to k z . 



the anisotropic s-wave pairing potentials. The graphical 
representations of these different pairing potentials are 
illustrated in Figs. |](a)-(d). 

Among different Ai g -rcprcscntations, the lowest-order 
possibilities include the isotropic s-wave order parame- 
ter, anisotropic s-wave with 6-fold in-plane modulations, 
or anisotropic s-wave with uniaxial symmetry, with the 
latter two illustrated in Figs. |l|(a) and[l](b). The lowest- 
order ^^-representation consists of twelve "lobes" of 
alternating phases, and the phases are even under k z - 
inversion. For either B\ g or _B2 S -representation, the or- 
der parameter consists of twelve lobes with alternating 



phases, and the phases are odd under fc 2 -inversion. 

To obtain the quasiparticle spectra for all possible 
pairing channels with different A&, we consider a crys- 
talline plane with a normal vector hk characterized by 
the parameters {9k,4>k)- Defining the direction of an 
incident quasiparticle relative to hk by the parameters 
(■{6im4>in) = (in, which explicitly considers the trans- 
verse momentum for the incident quasiparticles (i.e., a fi- 
nite "tunneling cone" ) relative to hk , such that 9i n is pri- 
marily confined between —(3 and /3, and < <pi n < 2tt, we 
can generalize the theory of Blonder-Tinkham-Klapwijk 
(BTK) (If, H to three dimensions (3D), and com- 
pute the tunneling current Ins as a function of the bias 
voltage V, temperature T, tunneling barrier strength Z, 
tunneling direction hk, and tunneling cone /3: 

Ins = Gnn / d(f) in \ d9 in cos 9 in e dE k 
Jo Jo J 

x[l+A-B]x [f(E k - eV) - f(E k )} . (2) 

In Eq. (0), Gnn denotes the normal-state conductance, 
Ek is the quasiparticle energy, A and B represent the 
kernals for Andreev and normal reflection, respectively, 
and f(E k ) is the Fermi function [ff, |§. Thus, the dif- 
ferential conductance spectra (dlNS /dV)-vs.-V can be 
obtained for given hk and A^ using Eqs. (0) and @. 
The representative spectra for high-impedance tunneling 
barrier Z = 5 are shown in the left panels of Fig. |l|(a)-(d). 

Except for the ^4i g -representation, the spectral char- 
acteristics for all other representations exhibit strong di- 
rectionality (i.e., dependence on the crystalline normal 
hk relative to the average quasiparticle momentum), as 
manifested by calculated spectra in the right panels of 
Fig. |j](a)-(d). It is clear that the ZBCP would have been 
a common occurrence in the tunneling spectra of MgB2 
pellets had the order parameter been one of the uncon- 
ventional pairing channels (A2 g , Bi g , B-2 g ). 

To compare the calculated results with experiments, 
we performed scanning tunneling spectroscopy on high- 
density pellets |2l], |H ^3| and c-axis textured films of 
MgB 2 f|§ at 4.2 K. Both the pellets and c-axis films 
were fully characterized |l2|, plj |2^, |24|, showing 
single-phased material with superconducting transition 
at T c = 39.0 K, sharp magnetization transition widths 
(AT C < 1 K for the pellets and AT C - 0.7 K for the 
films), and nearly 100% bulk superconducting volume 
lHI, ||, ||l |4}. According to XPS studies on these 
samples fl2|] , the surface MgC03 and Mg(OH)2 impu- 
rities on the as-grown MgB2 could be mostly removed 
by chemical etching, with no discernible etch residues for 
the tunneling experiments ]l2| . Tunneling studies were 
conducted on the as-grown and etched MgB2 pellets and 
films at 4.2 K, using a low-temperature scanning tun- 
neling microscope. Spatially resolved tunneling spectra 
were taken on over 100 randomly oriented grains of each 
sample. On each grain, the spectra were taken under 
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FIG. 2: (a) Spatially resolved tunneling spectra of a high- 
density MgB2 pellet. The main panel and the lower right inset 
illustrate data taken at locations 10 ~ 15 nm apart within one 
grain after and before chemical etching, respectively. The 
upper left inset shows representative spectra on the etched 
pellet with different junction resistance at 20 mV: 1) 108 Mfi, 
2) 179 MO, 3) 253 MQ. The work function for these spectra 
is typically 0.1 ~ 1 eV. (b) A series of tunneling spectra 
on an etched c-axis film (main panel), showing long-range 
spatial homogeneity in the spectral peak-to-peak energies and 
a large junction resistance ~ 330 MQ, (inset), (c) An image 
of the surface topography of the etched sample over an area 
(196nm x 60nm). The full scale for the height is 4.7 nm. 



the vacuum tunneling condition and on an area approxi- 
mately (200nm X 200nm) in size with nano-scale surface 
flatness. A large number of grains were studied on each 
sample to ensure sufficient statistical sampling of differ- 
ent hk in pellets. 

Representative tunneling spectra for a MgB2 pellet af- 
ter etching are shown in the main panel of Fig. |](a), 
and those for the same sample before etching are given 
in the lower right inset. We note significantly improved 
spectra after etching, with long-range spatial homogene- 
ity (> 400 nm) within each grain, which correlated well 
with the long-range atomic flatness of the topography 
as exemplified in Figure 2(c) and also according to our 
AFM images, and was in contrast to the strong spa- 
tial variations in both the spectra and topography of 
MgB2 powder p5|. Furthermore, the density of states 
(DOS) nearly vanished at Ep, with a normalized value 
[{dI NS /dV) v= o/(dI N s/dV)v=20mev] ~ 2%. While the 
tunneling spectra were homogeneous within each grain 
(with lateral dimension ~ a few fim [ p2[ ) , the gap values 
varied from grain to grain in the pellets, ranging from 
~ 5 to ~ 8 meV. On the other hand, tunneling spec- 
tra of etched c-axis oriented films were homogeneous ev- 
erywhere. Overall, no ZBCP was observed among over 
five hundred spectra taken on all samples. We therefore 
conclude that the pairing symmetry must be of the Ai g 
representation. 



To identify the correct pairing potential under the Ai g 
representation in Eq. ([l]) , we performed the BTK analysis 
for both anisotropic and isotropic s-wave pairing, as well 
as the isotropic BCS fitting to all spectra. The latter 
involved a disorder parameter T for an isotropic gap A 
with the density of states N(E) given by p(J : 



■A/"(E)=Re (-E - iT)/y/{E - iY) 2 - A 2 



oc dlNs/dV. 



For both BTK and BCS isotropic s-wave fitting, we 
notice several difficulties. First, the inclusion of the 
disorder-induced pair-breaking strength T alone cannot 
fully account for the spectral characteristics, particularly 
the line- width and line-shape of the peaks, as manifested 
in Figs. ||(a) and ||(b). Second, significant variations in 
the supposedly isotropic pairing potential must be in- 
voked to account for all data taken on the pellets. The 
variation was unlikely the result of bulk stoichiometric in- 
homogeneity because of the sharp superconducting tran- 
sition width (< 1 K) revealed in the magnetization mea- 
surement of our MgB2 pellet, fn other words, had the gap 
variation been the result of the grain-to-grain stoichio- 
metric variation, we would have observed a very broad 
T c distribution in the magnetization measurements, from 
~ 39 K to ~ 24 K for the 5^8 meV gap variation. 
Given the quality of the spectra and topography of our 
well characterized sample surfaces, we suggest that the 
variation observed in the gap values of MgBi pellets is 
the result of different grain orientations relative to the 
incident quasiparticles. The single gap value in the c- 
axis oriented films further corroborates the notion of k- 
dependent pairing potential. More importantly, had the 
pairing symmetry been isotropic s-wave, the (2A/fceT c ) 
ratios deduced from our tunneling spectra would not have 
varied from ~ 2.5 to ~ 4.5 from grain to grain for T c 
variation smaller than 1.0 K. In addition, to date there 
is no known theory for isotropic s-wave superconductors 
that can justify a (2A/fcsT c ) ratio smaller than the BCS 
value. 

On the other hand, the electronic and structural 
anisotropy in the MgB2 system can lead to anisotropic 
s-wave pairing, and therefore a fc-dependent pairing po- 
tential and a range of gap values in the STS studies 
of polycrystalline samples. Comparing the two possi- 
bilities of anisotropic s-wave pairing potentials depicted 
in Figs. |l](a) and |l|(b), we note that an in-plane 6- fold 
anisotropy would have resulted in a c-axis spectrum with 
a sharp peak at Aq(1 + e) and complicated spectral cur- 
vatures in M{E) for A (l — e) < E < A (l + e), as shown 
in Fig. |l](a). Such behavior was never seen in our data. 
In contrast, spectra derived from the order parameter in 
Fig. ^(b) appeared to be most consistent with our finding 
of smooth spectra on all samples, and with one maximum 
gap value at A z ~ 8 meV for the c-axis films. 

Using the anisotropic s-wave pairing potential 
Afe = A xy sin 2 9k + A z cos 2 9k, with the minimum gap 
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FIG. 3: BTK anisotropic and isotropic s-wave fitting, to- 
gether with the isotropic BCS fitting to representative spectra 
of (a) an etched MgE>2 pellet, and (b) an etched c-axis film. 
Given empirical values of A xy and A z , the anisotropic s-wave 
fitting is only sensitive to the variation in 9k and is insensitive 
to a wide range of (3 values that we have tested, from (71-/I8) 
to (7r/2). The fitting curves shown have assumed the most 
general case with f3 = tt/2. 

A xy « 5 meV and the maximum gap A 2 8 meV de- 
termined empirically, we can consistently account for all 
experimental data on both pellets and c-axis films by 
varying one parameter 9k- As exemplified in the main 
panel and inset of Fig. ||(a), the former is consistent with 
Ok = (tt/5) and the latter with 9k = 0. Similarly, the 
same pairing potential can also be applied to the c-axis 
film data with 9k = 0, as shown in the main panel of 
Fig. ||(b) . Our empirical finding of a smaller in-plane gap 
value {A xy < A z ) is consistent with the stronger in-plane 
Coulomb repulsion in MgB2 ^7j. A similar anisotropic 
s-wave pairing scenario has also been proposed recently 
to account for the thermodynamic and optical properties 
of MgB2 wires [ p8| . Furthermore, a number of recent 
experimental reports, including the upper critical field 
(H c2 ) measuremen ts |29[ |3(| , high-resolution photoemis- 
sion spectroscopy and electron spin resonance 
are supportive of significantly anisotropic properties in 
the superconducting state of MgB 2 . 

In summary, we have investigated the possible pair- 
ing channels in MgB2 based on group theory considera- 
tion, and have calculated the quasiparticle spectra using 
a generalized BTK theory for quasiparticle tunneling in 
3D. Comparing the calculated results with spectra taken 
on fully characterized MgB2 pellets and c-axis oriented 
films, we conclude that the order parameter of MgB2 be- 
longs to the Aig-representation of Dq^ group, and is best 
described by an anisotropic s-wave pairing potential with 



uniaxial symmetry. 
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